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Abstract —Magnetic induction (MI) based communication and 
power transfer systems have gained an increased attention in 
the recent years. Typical applications for these systems lie in the 
area of wireless charging, near-field communication, and wireless 
sensor networks. For an optimal system performance, the power 
efficiency needs to be maximized. Typically, this optimization 
refers to the impedance matching and tracking of the split- 
frequencies. However, an important role of magnitude and phase 
of the input signal has been mostly overlooked. Especially for 
the wireless power transfer systems with multiple transmitter 
coils, the optimization of the transmit signals can dramatically 
improve the power efficiency. In this work, we propose an 
iterative algorithm for the optimization of the transmit signals 
for a transmitter with three orthogonal coils and multiple single 
coil receivers. The proposed scheme significantly outperforms the 
traditional baseline algorithms in terms of power efficiency. 

I. Introduction 

Magnetic induction (MI) based transmissions are well 
known in the context of near-field communication (NFC) 
m, wireless power transfer (WPT) El, and wireless sensor 
networks (WSN) in challenging environments 0, 0. In this 
work, our main focus lies in the WPT using resonant coupling 
of magnetic antennas. Typically, MI based power transfer is 
only useful within a short range due to a dramatically low 
power efficiency otherwise, as it has been confirmed in nu¬ 
merous previous works (e.g. 0). Furthermore, the alignment 
of coils has a strong impact on the transfer efficiency, see 
e.g. 0, 0. Several attempts have been made to extend 
the magnetic induction based point-to-point transmission to a 
system with multiple receivers 0, multiple transmitters [9], or 
even multiple relays mol, ED- Moreover, MI based networks 
with multiple transceivers and relays have been analyzed 
for the underground WSNs 02. Finally, the multiple-input 
multiple-output (MIMO) technique has been introduced for 
different constellations of MI based communication and WPT 
systems, see e.g. El, 03. In particular, 03 suggests the 
use of a transmitter equipped with three orthogonally deployed 
coils. This enables a steerable directionality of the resulting 
magnetic field, the so-called magnetic vector modulation. 
This decomposition is based on the explicit weighting of 
the field vectors given by the orientations of the transmitter 
coils El- Unfortunately, this approach does not take into 


account the influence of multiple receivers on each other or on 
the transmitter. Similarly, in fl3l , only the carrier frequency 
is optimized, leaving the choice of e.g. signal phase and 
amplitude (beamforming coefficient^) suboptimal. 
Beamforming is a well known technique for maximizing the 
power efficiency of a MIMO system. In traditional radio 
frequency (RF) systems, the power efficiency optimization 
corresponds to the maximization of the receive power for 
fixed /. 2 -norm of the beamforming vector, because the con¬ 
sumed power in the transmitting device depends only on the 
constraint on L 2 -norm of the beamforming vector, not on 
its particular coefficients. For MI based WPT, the efficiency 
depends explicitly on the coupling between the transceivers, 
such that improving the coupling yields an increase of the 
power efficiency 0. Furthermore, since the signal reflection is 
approximately proportional to the squared mutual inductance 
El, an influence of the reflected signals on the transmit power 
is inevitable. In particular, these reflected signals can overlap 
constructively or destructively depending on the phase of the 
input signals. Hence, the transmit power depends on the choice 
of the beamforming coefficients and a sole receive power 
maximization becomes insufficient. Therefore, we propose an 
iterative algorithm, which takes into account all couplings 
between the coils and maximizes the WPT efficiency. 

For this work, we consider one transmitter with three ortho¬ 
gonal coils free of self-interference and multiple single antenna 
receivers randomly deployed in the near-held of the transmit¬ 
ter. For a more flexible system design, the power transfer links 
can be assigned different priorities, such that more power can 
be steered into the preferred direction. In this context, high 
efficiency gains can be observed compared to the mentioned 
baseline schemes. In addition, a WPT efficiency of up to 96% 
can be achieved even for more than three receivers. 

This paper is organized as follows. Section lU provides insight 
into the system model for WPT with three transmission coils 
and multiple receivers. Our system model allows also a priority 
based optimization, such that some receivers may get more 

*We follow the convention of WPT using electromagnetic waves 0) and 
adopt the term beamforming for the optimization of the transmit signal vector 
in spatial domain. 



Fig. 1. Example of a WPT system with 3D-coil based transmitter ( Tx ) 
and multiple (K = 3) single coil receivers {Rx\, Rx 2 , and Rx 3 ). Different 
energy levels motivate the priority aware efficiency optimization. 



Fig. 2. A single WPT link. Transmitter and receiver resonance circuits. 


power upon request. In Section [Till the efficiency optimization 
problem is formulated and different solutions are introduced. 
In Section IIV1 numerical results are provided, and Section [V] 
concludes the paper. 


II. System Model 


In this work, we utilize one transmitter equipped with 
three orthogonally deployed coils (3D-coil) and multiple ( K) 
receivers with one coil each, see Fig. Q] where K = 3. Every 
coil with inductivity L is considered as part of a resonant 
circuit, which includes also a capacitor with capacitance C and 
a resistor with resistance R (modeling the copper resistance 
of the coil). The capacitance C is selected to make the 
circuit resonant at the resonance frequency /o = - - ^ LC - 
However, the actual operating trcquenc\@ / for the WPT 
can be subject to optimization. Furthermore, each receiver 
circuit contains a real-valued load resistor Z k , which can 
be optimized individually for different receivers in order to 
minimize the power reflection at the receiver according to the 
previous work El, see Fig. U Hence, the inner impedance of 
a resonant circuit with index k can be given by 

Zk,k = j27r/L + j2 n f(j R’ (1) 

if circuit k belongs to the transmitter, or 


■Z'fc.fc — j27r/L + 


1 

l^Jc 


+ R + ^L,kt 


( 2 ) 


if circuit k belongs to a receiver. The induced voltage is 
related to the coupling between the coils, which is determined 
by the mutual inductance M. The knowledge of the mutual 
inductance is very important for the optimization of WPT 


2 For the WPT, only one frequency is utilized 0. 


systems |2j] and can be determined either by channel estima¬ 
tion or by distance estimation under the assumption that all 
other system parameters (coil dimensions, polarization, etc.) 
are known. In this work, we assume that the exact value of 
the mutual inductance between any pair of coils is available 
to the transmitter. With this, the frequency selectivity of the 
MI channel is entirely known, such that not only the WPT can 
be established, but possibly also the information transmission. 
For this, only the signal pulse shape needs to be adaptccQ. 
However, MI based information transmission is beyond the 
scope of this work. 

The orientations and alignments of the coupled coils have a 
strong impact on M and correspondingly on the path loss, cf. 
0. Hence, we model the mutual inductance between coils k 
and Z by 

Mk,i = M k ,i ■ Jk,h (3) 

J k ,i = 2 sin 9 k sin 9i + cos 6 k cos 9i cos cj>, (4) 

cf- 0, where 6 k and 9i are the angles between the radial 
directions of the coils k and Z, respectively, and the line 
connecting the two coil centers. <f> is the angle difference 
between the coils’ axes in the plane, which is orthogonal to 
the direction of transmission. M k .i represents the (absolute) 
value of the mutual inductance for the case J k j = 1. For the 
following, we define Z k j = j2nfM k j, Vfc ^ Z. 

We consider the complex-valued amplitudes U k and I k of 
the voltages u k (i) = U k ■ eand currents i k {t) = I k • 
Vfc, respectively. For each coil k, the current amplitude 
I k in the resonant circuit depends on the current amplitudes 
VZ ^ k in all surrounding circuits via the voltage equation 

h ■ Z k , k + (h ■ Z ktl ) = U k , (5) 

l^k 

where U k is the complex-valued amplitude of the input volt¬ 
age. In the following, the superscripts (-) T and (•) H denote 
transpose and Hermitian transpose, respectively. We assign the 
first three coil indices to the transmitter and the remaining 
indices to the receivers. In order to calculate the currents in 
all circuits of the coupled network, a set of voltage equations 
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needs to be solved. Here, Ur.r is the complex-valued input 
voltage vector at the transmitter. Furthermore, 0 stands for the 
all-zero vector and represents the input voltages u k (t) = 0 for 
all coils k belonging to a receiver, since only the transmitter 
is supposed to generate power. l' Tx and \ c Rx denote the current 
vectors of the transmitter and the receiver circuits, respectively. 
The matrices Z Rx , Zch and Z Rx contain complex impedances 
and are defined in the following. In this work, we use a 30- 
coil based transmitter, which means that all three transmitter 

Tor information transmission, bandlimited signals like root-raised cosine 
pulses are typically utilized. For WPT, the transmitted signal is usually a sine 


wave. 




















coils’ axes are orthogonal to each other, such that 


Z Tx = 


0 0 

0 -^2,2 0 

0 0 Z^,3 


(7) 


holds. The receiver coils are not necessarily orthogonal and 
we obtain 


Z r x = 


Zaa 


Zi,K+3 


Zk+ 3,4 
Z K+3,K+3. 


( 8 ) 


The purely imaginary matrices Z ch and Z ch in © stand for 
the influence of the receiver coils onto the transmitter coils 
and vice versa, respectively. Hence, Z ch is defined by 


^ 4,1 

Z$a 

Zk+ 3,1 

^4,2 

Zb ,2 

Zk+ 3,2 

Z 4.3 

Zb,3 ‘ ‘ 

Zk+ 3,3. 


By inverting the impedance matrix in © using (lSJ, we obtain 
similar to © 


I Tx = (Ztx - ZcfcZ^Z^) U Tx - AU Tx , (10) 
I k = -Z^Z^AUr* = CV Tx , (11) 


with implicit definitions of A and C. As known from the 
fundamentals of electric power generation and transmission 
(e.g. QT), ED), in order to produce enough active power in 
electric circuits, the transmitter/generator needs to release also 
the reactive power, which corresponds to the imaginary part 
of the generated complex power. The reactive power is not 
absorbed by the load, but fluctuates between the power source 
and the load impedance. Furthermore, without the reactive 
power, the induction coils cannot be operated. Hence, both 
real and imaginary parts need to be taken into account, such 
that the magnitude of the generated complex power (the so- 
called apparent power) is a better reference for the maximum 
transmit power than the pure active power fl8l . In the fcth 
transmitter circuit, the apparent power is given by 


Pt,k = \U k h\ = \Uk\ |4|, (12) 

where |-| denotes the element-wise absolute value operator. 
Therefore, we define the total power provided by the trans¬ 
mitter as 

3 

ft,total = P t ,k = \Utx\ T |AUt*| . (13) 

fc=1 

For the received active power at the load resistor Z^i of the 
receiver circuit l we obtain 


P r ,l = \h \ 2 Z L>1 . (14) 


In practical WPT systems (in particular related to the sensor 
networks), the energy consumption in receiver devices may 
be not the same, such that the depletion rate of the batteries 
varies. Therefore, some receiver devices may require more 
power than the others, see Fig. Q] Hence, we introduce the 


priority coefficients Wi, VZ. With this, the weighted received 
power is given by 


K+3 

Pr,total = ^ ^ IT/ I j J 
1=4 

W 4 Z LA 


0 


0 • • • Wk+3Zl,K+3. 

where (fill and (fl4l> have been used. 


CU Tl , (15) 


III. Beamforming 

In this section, the design of the optimal input vector U Tx 
(MI beamforming vector) for the weighted power efficiency 
maximization is discussed. First, the optimization problem is 
formulated. Then, some of the most promising approaches 
including the proposed iterative algorithm are presented. 


A. Problem formulation 

In this work, we define the WPT efficiency with respect to 
the sum apparent power given by dT3l) . Hence, the optimization 
problem can be formulated as 


P, 


max 


r,total 


u T*,f,Z L ',Vl P(, total 


(16) 


where the parameters / and Zlj, V/ can be optimized ac¬ 
cording to the literature 0, ED. Hence, we focus on the 
optimization of U Tx- 

It can be shown that the problem (fl6l > is non-convex due to 
the non-convexity of Pt,total. suc h that the well-known convex 
optimization tools [ fl9l cannot be used. Therefore, we provide 
some suboptimal schemes. Our proposed algorithm discussed 
in Section IIII-B3I can be shown to reach a local optimum in 
case of convergence. 


B. Proposed algorithms 

In the following, three different approaches are described. 

1) Closest neighbor based beamforming: As known from 
the near-field communication and coupled-mode theory, usu¬ 
ally only the closest neighbors tend to establish a strong 
coupling, such that the path losses between the transmitter 
and any other receiver are dramatically larger, especially in 
case of weak couplings between coils. Based on this principle, 
the idea of optimizing the beamforming vector for the closest 
neighbor of the transmitter is motivated. Hence, using ©-©, 
and assuming M/ j 1 = Mip, = the magnetic induction 

at the receiver l is related to 

Ui tx [Ii,l2,13 ] J/,2, Ji, 3 \ T ‘ (17) 

where the current contributions from other receiver coils are 
neglected due to a very high path loss of such signal^ 
Furthermore, due to identical sets of circuit elements in all 
transmitter resonance circuits, Ik oc Ut■ k = {1,2,3} holds. 


4 Since here the target receiver l is the closest receiver to the transmitter, 
the signals arriving from the neighboring coils are basically heavily attenuated 
reflections from the receivers that are further away from the transmitter. 













which results directly from {7} and ( flOl ) for weakly coupled 
coils. Hence, we substitute this result in O and obtain 

© cx [UuU^UaliJi,!, ■ (18) 

Due to (fl8l) . the beamforming vector 

Utx = Ji,2, Ji, 3 ] T ■ (IV) (19) 


3) Proposed iterative algorithm: The proposed idea is to 
approximate the total transmit power by a squared L 2 -norm 
in each iteration of the algorithm. Using this approximation, 
the optimal beamforming vector is calculated, which helps up¬ 
dating the solution in the next iteration. For the approximation, 
we assume that in case of convergence of this algorithm, 

IUTx.nl « |U T x,n—1| (24) 


corresponds to the maximum ratio combining (MRC) solution. 
Obviously, this approach provides a close-to-optimum solution 
in case of a weak coupling between any adjacent coils and with 
a dominant receiver. Moreover, any variations of the matrices 
A and C due to a stronger coupling between neighboring 
coils result in a deviation of the optimum solution from this 
beamforming solution, thus decreasing the power efficiency. In 
addition, the transmit power is not taken into account at all, 
such that the power efficiency of this solution is suboptimal. 

2) Receive power maximization: The second approach is 
based on the eigenvalue decomposition of the receiver matrix 


D = C H 


WiZL,i 


0 


0 


c 


Wk+3^L,K+3_ 


( 20 ) 


and corresponds to a typical beamforming solution in the 
traditional RF systems. The eigenvector pertaining to the 
maximum eigenvalue maximizes the total receive power (ITD . 
This approach takes into account the impact of the coupling 
of coils onto the receive power. Hence, it provides a more 
accurate solution for the receive power maximization. 

For weak couplings between coils (low mutual inductance) 
and identical inner impedances Z ± t i = Zy> = Z 3 3 in the 
transmitter coils, matrix A is approximatel)lj given by 

A«Z^ = |Z lil r 1 -I, (21) 


where I is the identity matrix. The correctness of (ITH is 
confirmed by previous works, e.g. GH- This yields 

p t , total « U^Utx • |Zi,ir\ (22) 


where |Utx| T |Utx| = U^Urx h as been used in (ITD . 
Therefore, the beamforming optimization problem reduces to 
an eigenvalue problem given by 


max 

U Tx 


Ufxg Utx 

U^Utx 


(23) 


Hence, this approach is optimal for the weak couplings be¬ 
tween coils. With increasing mutual inductance, the approxi¬ 
mations (ED and ED are not valid anymore. Correspondingly, 
the efficiency ED becomes suboptimal, because the transmit 
power according to (ITD is not explicitly considered. Thus, a 
more powerful algorithm is proposed in the following. 


holds, where U Tx,n denotes the state of the vector Utx at the 
end of the nth iteration. At first, we approximate the transmit 
power (ITD by 

Pt ,total = |UTx,n| T |AUrx,n| ~ |UTx,n-l| T |AUrx,n| > (25) 

such that the order of the transmit power with respect to U Tx,n 
reduce^. Then, we express |Urx,n-i| T as 

|Ur*,„-if= [1,1,1] V n , (26) 


using matrix V„ defined by 


V 


n 


)Ui,n-i\ 0 0 

0 \U 2 ,n-l\ 0 

0 0 |t/ 3 ,n-l| 


By inserting ED into ED and using (ED . we obtain 


(27) 


Pt,total « [1,1,1] v„ |AUTx,n| = [1,1,1] |V„AUrx,n| ■ 

(28) 

Moreover, ESI) can be transformed into a squared L 2 - 
norm. For this, we define S„ = V„A and approximate 
|Y„AU T x,n| = |S„U T x,n| from d28j by 

ISnUrx.nl ~ |S„Urx,n|®(|S„Urx,n| 0 |S„Urx,n-l|) , (29) 


where 0 and 0 represent element-wise vector multiplica¬ 
tion and division, respectively. Hence, by reformulating ED . 
|S„Urx,n| can be expressed as 


|S„Urx 


S„Urx 


0 


SnUrx.n— 1 


© 


(30) 


(2) 

where (-) w denotes element-wise square operator. By multi¬ 
plying (ED with a vector [ 1 , 1 , 1 ], the transmit power can be 
expressed as a squared L 2 -norm using ESI i: 


|Urx,n| T |AUrx.n| 


[SnUrx.n 0 \J |SnUrx,n—1| 


2 

2" 


(31) 


For the clarity of exposition, we denote |S„Urx,n-i| by vector 
G„ = [Gi,G 2 ,G 3 ] T . Using G„, the element-wise division in 
(ED can be formulated as a multiplication with a matrix Q ra , 
where 


Q 


VUT 1 0 0 

0 v^" 1 0 

0 0 y/UZ 


Then, we obtain 

SnUrx.n 0 (SnUrx.n—11 = QnSnUrx.n- 


(32) 


(33) 


5 For 


2tt/ ma x{M k i} 
k.l 


< R, Zch Z lll Z Ch in m is negligible. 


6 The order with respect to the complex-valued variable n is larger for 
|U T x,n| |AU Ta ., n | than for |AU Ta ,n|. 














Finally, by inserting ( l33l ) into ( [3Tb . we obtain 


\^Tx,n\ |AUTa;,ri 


|Q„S n U T *, 112 


n|| 2 - 


(34) 


Based on < l34b . the beamforming problem (fl6] > can be reduced 
to a generalized eigenvalue problem in each iteration 


U Tx,n = argmax 


Uf x D V Tx 


TJt 


U^(Q„S„) ff (Q„S n )U T: 


(35) 


where D is given by (l20l >. The solution to this problem is 
typically computed using a substitution 


X„ = ( Q n S n )\J Tx , (36) 

Vtx = (Q„S„)” 1 X„. (37) 



a) F=0.1 b) F=1 



c) F=10 d) F=100 


An eigenvalue decomposition is applied to the matrix 
((QnS„) -1 ^ D ^(Q„S n ) -1 ^ and the eigenvector with the 
maximum eigenvalue is picked as the optimal solution for X„. 
Then, using (l37l i. the optimal beamforming vector U r x ,n is 
calculated. This vector replaces Ftx.u-i in the next iteration. 
For the starting point Utx,o, we choose the eigenvector 
according to Section HII-B2I Hence, we start with the receive 
power maximization without taking into account the influence 
of matrix (Q„S„) and then improve the energy efficiency using 
the proposed approach. For a stopping condition, a maximum 
number of iterations or a minimum efficiency gain of the 
current iteration over the previous iteration can be used. 

The major benefit of the proposed algorithm is due to the 
adaptation of the beamforming to the changes of the transmit 
power. In case of convergence, the approximations (l24l >. < [29b . 
and correspondingly ( l34t are valid. Then, the maximization 
problem in ( IT6b becomes concave, and the equivalent min¬ 
imization problem becomes convex OH and can be solved 
using the generalized eigenvalue decomposition. Thus, the 
algorithm leads to a locally optimal solution in case of con¬ 
vergence. Unfortunately, the existence or absence of any other 
locally optimum solutions cannot be shown mathematically, 
such that the obtained solution is not necessarily globally 
optimal. However, as shown in Section [IV] the proposed 
solution performs well and reaches high power efficiencies. 

IV. Numerical Results 

In this section, we present numerical results for the MI 
beamforming optimization. In order to follow the convention 
of the WPT community 0, we consider a factor f)./ = 
'2ivj M k ,i , w hi c h corresponds to the product of the quality factor 

2 ^ L and the coupling coefficient ^44 between coils k and 
l, respectively, in our performance investigations. We assume 
that all receivers are placed at the same distance d from the 
transmitter, such that the coupling coefficient is identical for 
all transmitter-receiver links, which means Mk,i = M, k £ 
{1,2,3}, i € {4,, K + 3} and F k ,i = F, k e {1, 2, 3}, l £ 
{4,..., K + 3}. Since the distances between the receivers may 
vary, the mutual inductance between them differs from M. 


Fig. 3. Beamforming pattern for F = (0.1,1,10,100}. 

Assuming that the distance between two adjacent receivers 
with indices l\ and I 2 is di lt i 2 , Mi lt i 2 can be expressed as 

M W2 =m(^-) , (38) 

because the mutual inductance scales with the third power 
of the transmission distance m. For simplicity, the carrier 
frequency is selected to be equal to the resonance frequency, 
/ = foil In order to provide insight into practically relevant 
values for the factor F, we calculate it for the following 
scenario. Assume, a single turn rectangular air core coil of 
4 cm x 6 cm cross-section area made from copper wire of 
3 mm thickness is operated at fo = 125 MHz. Then, for a 
transmission distance d = 0.4 m, we obtain F ss 15 li20ll . 

We start with the visualization of the beamforming for a single 
receiver. For this, we assume that the receiver with randomly 
rotated coil is at first placed at 0° of the angular space of 
the transmitter and the optimal transmit signal \Jt x , optimized is 
found using the proposed algorithm. This signal is then used in 
order to evaluate the beamforming effect of this particular con¬ 
stellation. For this, the receiver is moved around the transmitter 
and its axis orientation is randomly rotated. For each point in 
the angular space we determine the power efficiency using 
Utx, optimized- This calculation is repeated for 1000 different 
constellations. The mean value of the resulting efficiency 
pattern is shown in Fig. [3] For comparison, we also show 
the results for a uniform beamforming vector Utx, uniform = 
[1,1,1] T • (1 V). This baseline scheme is motivated by fl4l . 
where it is suggested to maximize the field strength in the 
preferred direction by simply increasing the transmit power in 
the circuit with corresponding coil orientation. On the contrary, 
using Uxa:, u n if orm , a quasi-omnidirectional field propagation is 
supposed to be achieved. As we can see from Figs. [3}t) and 
b), this field propagation is on average not omnidirectional and 
provides a better efficiency for the receivers deployed at 45° 
and 225°, respectively. The power efficiency of the proposed 

7 In principle, the optimal carrier frequency according to CCD could be 
selected. 






























Fig. 4. Example on power efficiencies versus F of two randomly deployed 
receivers using different priority metrics. 

solution becomes more and more directional with increasing 
F and converges to a single peak. Hence, the power is steered 
only into the direction of the dedicated user, which is very 
beneficial, since almost no power is lost in the non-preferred 
directions and very limited interference is imposed on the other 
MI based communication systems. 

In order to visualize how the priority aware efficiency maxi¬ 
mization affects the receive powers, we show an example on 
the power efficiencies of two receivers for different priority 
metrics, see Fig. 0] These receivers are located at randomly 
selected positions on the circle in distance d from the trans¬ 
mitter and their coils are randomly rotated. For comparison, 
also the results of the closest neighbor based optimization are 
presented, where always the first receiver (Rx 1) is selected 
as the closest neighbor due to the equal distance of both 
receivers to the transmitter. Obviously, for equal priority of 
the receivers, the second receiver (Rx 2) obtains more power 
than the first receiver using the proposed solution. The sum of 
their efficiencies is obviously also larger than that of the closest 
neighbor based approach. By prioritizing the first receiver 
over the second (2:1 priority), the receive power at Rx 1 
becomes significantly larger. Interestingly, this efficiency even 
is superior to that of the closest neighbor based optimization, 
which is supposed to maximize the receive power for Rx 1. 
This is due to the couplings between coils and the transmit 
power that have been not taken into account in this approach. 
For this particular example, the priority increase by factor 2 is 
already sufficient in order to steer the power into the direction 
of Rx 1. In general, larger priority factors (10:1 or even 20:1) 
may be needed in order to obtain the desired performance of 
the beamforming. 

Finally, we show the mean achievable power efficiency for 
WPT to multiple receivers using the algorithms described in 
this work. For this, we consider five receivers and set all 
priorities to 1, such that no scaling of the receive powers is 



Fig. 5. Average power efficiency versus F of WPT to five randomly 
distributed receivers using different algorithms. 

performed. The mean power efficiency is then calculated by 
averaging over the efficiencies from 1000 different constella¬ 
tions, where each constellation corresponds to the realization 
of a random placement of the given number of receivers in 
distance d around the transmitter and a random orientation 
of their coils. Also, for each constellation, the orientation 
of the transmitter is randomly selected. The results for the 
average power efficiency versus F depicted in Fig. [5] show 
significant efficiency improvements of up to 37% for the pro¬ 
posed algorithm compared to the other solutions. As discussed 
in Section IIII-B1I the closest neighbor based approach does 
not take into account the couplings between receivers and 
the transmit power variations, such that it mostly performs 
worse or equal to the maximum eigenvalue based approach. 
As mentioned in Section IIII-B21 the maximum eigenvalue 
based beamforming is optimal for weak coupling, therefore its 
efficiency curve overlaps with that of the proposed solution for 
F < 1. For stronger couplings (F > 1), the proposed solution 
shows a much steeper increase of the efficiency, yielding large 
efficiency gains. The mean power efficiency reaches values of 
0.96 for large F, which corresponds to a power loss of only 
4%. This makes the proposed algorithm very promising. 

V. Conclusion 

In this work, a novel beamforming solution for the efficiency 
maximization problem in MI based WPT systems is presented. 
Due to the non-convexity of the transmit power metric, the op¬ 
timum solution cannot be determined analytically. Therefore, 
we presented three suboptimal approaches and discussed their 
differences. The third approach corresponds to our proposed 
solution. It utilizes several assumptions, that reduce the order 
of the optimization problem. In the convergence point of 
the algorithm, a locally optimal solution for this problem 
is obtained. Furthermore, a priority aware optimization is 
possible, that allows for allocation of the power to different 






















receivers according to their needs. For the WPT to multiple 

receivers, significant gains have been observed, which makes 

the use of our algorithm very promising. 
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